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Abstract 

Alpha elastic scattering angular distributions of the Q:)^°®Cd reaction were measured at three energies around 

the Coulomb barrier to provide a sensitive test for the a + nucleus optical potential parameter sets. Furthermore, 
the new high precision angular distributions, together with the data available from the literature were used to study 
the energy dependence of the locally optimized a + nucleus optical potential in a wide energy region ranging from 
P,Lab = 27.0 MeV down to 16.1 MeV. 

The potentials under study are a basic prerequisite for the prediction of a-induced reaction cross sections and 
thus, for the calculation of stellar reaction rates used for the astrophysical 7 process. Therefore, statistical model 
predictions using as input the optical potentials discussed in the present work are compared to the available ^'^®Cd 
+ alpha cross section data. 

Keywords: astrophysical 7 process; elastic alpha scattering experiments, a + nucleus optical potential 


1. Introduction 

The bulk of the naturally occurring nuclei heavier than iron observed in the Solar System 
were synthesized via neutron capture processes. In the case of the so-called astrophysical 
s process, the neutron flux is moderate, elements are synthesized through slow neutron 
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captures and (3 decays up to 0,3. On the contrary, in an astrophysical environment 
where the neutron densities are orders of magnitude higher than the ones available for the s 
process, the neutron captures become faster than the /d-decays and extremely neutron-rich 
nuclides close to the neutron drip line can be synthesized. These nuclei decay back into the 
valley of stability when the neutron flux ceases. This process — the so-called astrophysical 
r process — can only take place in explosive environments 0 , 0 . 

Furthermore, there are about 30-35 proton-rich nuclides between ^^Se and ^®®Hg, which 
cannot be formed by neutron capture processes 0,33- They are the so-called p isotopes. 
The number of p isotopes depends on the one hand on the state of art s process models: a 
recent investigation 0] showed that there are large s process contributions to the abundance 
of ^^'^Er, ^^^Gd and ^®°Ta. On the other hand, the precise estimation of the r process 
contribution to the abundance of certain p isotopes (e.g. ^^^In or ^^^Sn, see in 0 ) is essential, 
too. Consequently the modihcation of the list of p isotopes could be necessary. 

The natural isotopic abundance of the p isotopes is typically 10 to 100 times less than 
that of the more neutron-rich isotopes of the same element that were created in the s or r 
processes. In their production, photon-induced reactions at temperatures around a few GK 
play a crucial role. It is generally accepted that the main stellar mechanism synthesizing 
these nuclei — the so-called 7 process |l 0 | — is initiated by photodisintegration reactions 
on pre-existing neutron-rich s and r seed nuclei. Photons with high energy and high flux, 
necessary for the photodisintegrations, are available in explosive nucleosynthesis scenarios, 
as in the Ne/0 burning layer in Type II- 0,33 or in the Type la supernovae jll|, flil ]. 
However, based solely on this model, the calculations are not able to reproduce the Solar 
System p abundances mainly in two mass regions: 

— The light p isotopes near the Mo-Ru region are largely underproduced by the pho¬ 
todisintegration reactions. In recent years, several other processes have been discussed that 
could contribute to the production of light p nuclei through proton captures in a hi ghly 


proton-rich material: the so-called rp-(rapid proton capture) process [l3|, i/p-process [l4 


and proton c^tures on highly enriched s process seeds during Type la supernova explo¬ 


sion 


111 . I 12 I. Il5| . However, the role of these processes in nucleosynthesis remains unclear. 


e.g., for the rp-process it is uncertain whether the produced nuclei can be expelled into the 
interstellar medium and plays any role in the galactic chemical evolution. 

T he p nuclei in the intermediate region between 150 < A < 165 are underproduced. 


too [16|, II 2 I. It remains unclear whether this dehciency is due to nuclear cross sections, 
stellar physics, or if alternative / additional process has to be invoked. 

Although the basics of the 7 process were laid down several decades ago, many details 
of the process are still unknown. On the one hand, this concerns the ambiguities in the as¬ 
trophysical conditions under which the process takes place (seed isotope abundances, peak 
temperatures, time scale, etc). On the other hand, large uncertainties are introduced into 
the calculations by the nuclear physics input, most importantly by the reaction rates (deter¬ 
mined from cross sections). All 7 process models require the use of huge reaction networks 
including tens of thousands of nuclear reactions, and the rates of these reactions at a given 
stellar temperature are necessary inputs to the network calculations. The reaction rates are 
generally taken from calculations using the Hauser-Feshbach (H-F) statistical model 

2 


18, 19 




















/ Nuclear Physics A 00 (2015) iPT^ 


3 


The accuracy of the H-F predictions mainly depends on the adopted nuclear models for 
optical-model potentials, y-ray strength functions, and nuclear-level densities. Previous 7 
process related experiments showed that above the A 100 region the precise knowledge 


on the a -|- nucleus optical potential would be highly desirable 
order to increase the predictive power of the model calculations. 


20|, I21I, l22|, l23|, 124, l25 


m 


The subject of the present work is the study of the angular distributions of elastically 
scattered alpha particles on the ^°®Cd p nucleus. On the one hand experimental angular 
distributions measured at higher energies are available in the literature and therefore, 
the opportunity to study the energy dependence of the imaginary part of the a -|- nucleus 
optical potential in a wide energy region ranging from 'Ei^ah = 27 MeV down to ^Lah = 16.1 
MeV is given. Furthermore, recently the cross sections of the a-induced reactions on ^°®Cd 
were measured at very low energies, close to the astrophysically relevant ones 


27 . It was 


found that there is about a factor of 3 difference between the experimental cross sections 
and the H-F predictions in the case of the (a, 7) reaction. The inaccurate knowledge on the 
a + nucleus optical potential was identihed as the source of this disagreement. 

Finally, the elastic scattering data obtained extends the already existing data set of 
experiments performed on several proton-rich nuclei 26, 28, 29, 30, 31, 32, 33, 34, 35, 36| and 
provides additional data to test the global parametrizations used in astrophysical network 
calculations. 

The paper is organized as follows. Experimental details are described in Sec. [H The 
measured angular distributions are presented in Sec. El Here, the new data are analyzed 
using a locally adjusted potential, and predictions from global a-nucleus potentials are 
shown. In Sec. 01 we compare the results of the H-F calculations — calculated with the 
use of the local and global potentials discussed in the present work — with the available 
a-induced cross section data taken from the literature 
in Sec. 01 


27|. Finally, conclusions are drawn 


2. Experimental procedure 

The experiments were performed at the cyclotron laboratory of Atomki, in Debrecen, 
Hungary. The alpha elastic scattering angular distributions between 'diab = 20° and 'Oiab = 
170° were measured on the p nucleus ^°®Cd at three energies around the Coulomb barrier. 
Highly enriched ^°®Cd targets (~ 97%) of 250 pg/cm^ thickness were prepared using vac¬ 
uum evaporation technique. The experimental setup was similar to that used in previous 


experiments 36|] (and references therein). In the following paragraphs, a short outline of 


the most relevant characteristics of the experimental setup is given. A detailed description 
of the experimental setup can be found in e.g. (j^ . 

The angular distributions of the elastically scattered alpha particles were measured with 
the help of two pairs of silicon detectors, mounted on two independently rotatable turntables. 
The solid angles of the detectors were between 1.63 x 10“^ sr and 1.55 x 10“^ sr. In addition, 
two silicon detectors under 'diab = 15°, with a solid angle of 8.10 x 10“® sr were used for 
absolute normalization, to monitor the target thickness and to take into account the beam 
position on the target. 
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The energies of the alpha particle beam were E^ab = 16.1 MeV, 17.7 MeV and 19.6 MeV 
with typical beam intensities in the order of 250 nA. Two different sets of spectra registered 
during the measurement of the elastic scattering cross section are shown in FiglU On the 
left side of the figure, the spectra measured at position Oiab = 25° are shown for the three 
measured energies. The right side of the figure presents the spectra recorded at E^ab = 19-6 
MeV, for three different angular positions. As it can be seen, the elastic peak can be easily 
separated from the additional components, both at forward (carbon and oxygen peaks) and 
backward (inelastic scattering on ^°®Cd) angles. In addition, pulser signals were fed into 
the pre-amplihers in order to have a control on the dead-time of the detection system. The 
pulser signal amplitude was set so that the pulser peak always falls above the ^°®Cd elastic 
peak (see Fig. [T]). 

Knowledge of the exact angular position of the detectors is of crucial importance for the 
precision of a scattering experiment because the Rutherford cross section depends sensitively 
on the angle. The uncertainty in the Rutherford normalized cross section is dominated 
by the error of the scattering angles in the forward region. In a similar way to previous 
experiments 36| (and references therein) we took advantage of the kinematics of the reaction 
^^C(a, a)^^C and measured the elastically scattered alpha particles in coincidence with the 
recoil We hxed one of the detectors on one side of the chamber at i9a = 80° to measure 
the alpha particles, and measured in coincidence the recoil nuclei around the expected 
angle of 1 ) 12 ^ = 40.38°. The result can be seen in Fig. [2J The uncertainty of the scattering 
angle was determined to be Ai9 = 0.07°. Several measurements were performed at other 
angular positions, and confirmed this result, ensuring the high precision of the experimental 
data. 


3. Results and discussion 

Complete angular distributions between 20° and 170° were measured at energies of E^ab 
= 16.1 MeV, 17.7 MeV and 19.6 MeV in 1° (20° < i) < 100°), 1.5° (100° < < 140°), 

and 2° (140° < < 170°) steps. The statistical uncertainties varied between 0.3% (forward 

angles) and 1 to 2% (backward angles) j^. The count rates V('d) have been normalized to 
the yield of the monitor detectors NuonX’^ = 15°): 

f (a) - ( 

[dnJ^^^NMon An ^ ’ 

The solid angle of the detectors is represented by All. The Rutherford-normalized cross 
sections cover around 2 orders of magnitude, taking into account all the angular distributions 
measured {Eiab = 16.1 MeV to ELab = 27.0 MeV). Over this range of cross sections, our 
measured angular distributions have an accuracy of around 4 to 5%, while the data in ( 2 ^ 
claim around 1 to 2% at forward angles and the authors add < 10% for backward angles 
due to the relative low enrichment. 

In the following, the analysis of the new experimental data in the framework of the 
optical model will be presented. Our analysis is extended up to E^ab = 27 MeV by taking 

into account the elastic scattering angular distributions from ( 2 ^ . 
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Figure 1. (Colour online) Typical spectra taken for the reaction ^°®Cd(Q;, a)^°®Cd. On the left panels, the 
spectra taken at a fixed angle drab = 25° are displayed for all three measured energies. The right side of 
the figure shows spectra at forward, middle, and backward angles for a fixed energy pLab = 19-6 MeV. The 
elastic peak can be clearly distinguished from other peaks in all cases. In each spectrum, the pulser signals 
used to control the dead time were set at energies higher than the elastic peak. The most prominent peaks 
are labeled with the corresponding excitation energies of the residual nuclei. 


The optical model potential combines a Coulomb term with the complex form of the 
nuclear potential, which consists of a real and an imaginary part. Usually, the parameters of 
the optical potential are derived from the analysis of the angular distributions of elastically 
scattered alpha particles and are adjusted to experimental a-induced cross sections if they 
are known. The a + nucleus optical potential is given by: 
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Angular Calibration 



Figure 2. Results of the coincidence measurement. The relative number of recoils measured in 
coincidence with the elastically scattered alpha particles (da = 80°) is shown as a function of the recoil 
angle. The gray region shows the expected position for the maximum rate. The solid line is a Gaussian 
curve to guide the eye. The accuracy of the measured angle was determined to be Ar? = 0.07°. 


U{r) = Vc{r) + V{r)+iW{r) (2) 

where Vc(r) is the Coulomb potential (calculated from a homogeneously charged sphere), 
V (r) and W (r) are the real and imaginary parts of the nuclear potential, respectively. In the 
following, we will hrst determine a locally adjusted potential (Sec. 13.ip and then compare 
the new experimental data to predictions of global a + nucleus optical potentials fSec. 13.2p . 


3.1. Local optical potential analysis 

A local optical potential is tailored to describe the experimental data of a small region of 
the mass range, i.e., it is a potential with parameters derived from the data of a few alpha 
elastic angular distributions or even from just a single experiment. This type of potential is 
used as an analysis tool, since it is possible to achieve Xred close to unity, and therefore 
it is possible to study the impact/sensitivity of its different parameters in the description of 
the experimental data. 

The real part of the local potential is based on the double folding procedure of js^ 
calculated using charge densities based on electron scattering js^; this V{r) potential is 
modihed by a strength parameter A and a width parameter w, as in the case of the global 
potential ATOMKI-Vl (see Sec. 13.2p . The nucleon-nucleon interaction is parametrized using 


the well known DDM3Y interaction |23,1^, l4]J. The htted local potential has a Woods- 
Saxon form imaginary part lT(r), which consists of a surface Ws term only at very low 
energies 
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Table 1. Parameters of the local optical potentials derived from the fits to the elastic scattering angular 
distributions measured at Atomki (Atomki) and at University of Notre Dame [1^ (ND): for further details, 
see text. 


Ec.m 

A 

W 

Jr 

^rms,R 

Ws 

Rs 

as 

Ji 

^rms,I 

^reac 

Xied 

(MeV) 

(-) 

(-) 

(MeV fm^) 

(fm) 

(MeV) 

(fm) 

(fm) 

(MeV fm^) 

(fm) 

(mb) 

(-) 

15.6 - Atomki 

1.218 

1.029 

351.84 

5.3602 

193.01 

1.563 

0.291 

91.66 

7.490 

341.9 

0.8 

16.4 - ND 

1.211 

1.020 

341.06 

5.3152 

169.47 

1.468 

0.413 

101.19 

7.1421 

491.3 

1.0 

17.0 - Atomki 

1.373 

1.004 

368.91 

5.2324 

166.30 

1.464 

0.403 

96.51 

7.1158 

544.4 

0.9 

17.0 - ND 

1.423 

1.007 

385.38 

5.2462 

171.17 

1.446 

0.436 

105.02 

7.0679 

581.2 

0.7 

18.3 - ND 

1.371 

1.005 

369.37 

5.2370 

149.35 

1.428 

0.464 

95.30 

7.0147 

733.4 

1.1 

18.9 - Atomki 

1.378 

1.000 

365.29 

5.2076 

138.73 

1.402 

0.505 

93.26 

6.9422 

803.0 

0.8 

21.2 - ND 

1.380 

1.008 

374.43 

5.2500 

134.37 

1.419 

0.500 

91.41 

7.0108 

1039.5 

2.0 

21.2 - ND* 

1.393 

1.004 

373.74 

5.2304 

137.40 

1.414 

0.492 

91.26 

6.9812 

1020.9 

1.7 

Average 

1.336 

1.010 

365.18 

5.2641 

160.34 

1.456 

0.430 

96.34 

7.1119 

- 

- 

Av. w/ Ren. 

1.338 

1.010 

365.08 

5.2613 

160.78 

1.455 

0.429 

96.31 

7.1077 

- 

- 

23.6 - ND 

1.401 

0.994 

364.39 

5.1760 

159.16 

1.372 

0.492 

99.72 

6.7903 

1150.6 

5.4 

23.6 - ND* 

1.354 

1.007 

366.90 

5.2475 

152.59 

1.388 

0.513 

102.18 

6.8891 

1209.2 

2.6 

26.0 - ND 

1.446 

0.989 

370.83 

5.1515 

132.06 

1.372 

0.512 

86.21 

6.8157 

1287.2 

13.2 

26.0 - ND* 

1.392 

1.004 

373.17 

5.2288 

128.98 

1.392 

0.525 

88.96 

6.9209 

1346.5 

7.4 

* - Abs. renormalization; 1.016 at 21.2 MeV ; 0.946 at 23.6 MeV ; 0.929 at 26.0 MeV 

23.6 - S+V t 

1.433 

0.990 

368.14 

5.1547 

7.65 

1.686 

0.423 

54.10 

4.8021 

1153.2 

1.4 


= 

31.45 



Rv — 

1.131 



= 0.145 


26.0 - S+V t 

1.413 

0.996 

370.10 

5.1880 

10.65 

1.648 

0.464 

54.82 

5.0561 

1309.1 

2.5 


= 

29.98 



Rv = 

1.144 



a-v = 0.162 



21.2 - Phase Shift 1006.4 1.2 

23.6 - Phase Shift f 1154.7 1.1 

26.0 - Phase Shift f 1309.1 1.2 


t - Abs. renormalization: 0.959 at 23.6 MeV ; 0.921 at 26.0 MeV 


By fitting the data from each of the angular distributions, the Ws, Rs and as parameters 
of the imaginary part were left as free parameters for the htting. For the real part the bare 
double folding potential was calculated by the DFOLD code , using the aforementioned 
DDM3Y interaction, and was hxed for all the Eiab- The width w and strength A parameters 
of the real part were left free, in order to modify the bare folding potential (with A = tc = 1). 

The results of this local potential study are shown in Table [U The table lists all the 
resulting parameters, as well as the calculations for the volume integrals Jr and Jj, root 
mean-square radii rrms,R and Vrmsj, the reaction total cross section areac and the obtained 
Xled- 


It is obvious from Table [T] that the angular distributions above 20 MeV measured at 
Notre Dame can only be described with relatively poor x^ed- There are two reasons for 
this unexpected large 

First, the absolute values of the Notre Dame data have been normalized to Ruther¬ 
ford scattering {(t/ctr = 1.0) at the most forward experimental angle. However, this proce¬ 
dure is not justihed at the higher energies. We have allowed to vary the absolute normaliza¬ 
tion of these angular distributions, leading to Xred which are about a factor of two smaller 
than the ones using the original normalization. The normalization factors are 1.016, 0.946, 
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and 0.929 at 21.2 MeV, 23.6 MeV, and 26.0 MeV, respectively. 

Second, as already mentioned above, a surface imaginary potential is not sufficient at 
higher energies, and consequently a volume Woods-Saxon potential has been added (labelled 
“S+V” in Table [T]). This improves the obtained Xred ^y another factor of two. 

As a further test, in addition to the above mentioned local potential hts, phase shift hts 


were also performed for the Notre Dame data j26| using the method of j^. Here again 
Xred ~ 1 be obtained with modihed absolute normalization of the data, and the 

best-£t values for the absolute normalization agree within about 1 % with the local potential 
hts. This confirms that the absolute normalization of these data has to be revised. 

For the prediction of low-energy a-induced cross sections we have chosen an average po¬ 
tential from all angular distributions below 21.2 MeV which can be described using a pure 
surface Woods-Saxon imaginary potential. These average values are practically insensitive to 
the modified normalization of the Notre Dame data because the strongly affected data at 
23.6 MeV and 26.0 MeV do not contribute here (the additional volume Woods-Saxon imag¬ 
inary potential complicates the extrapolation to low energies). The minor renormalization 
of the 21.2 MeV data affects the average local potential parameters by the order of 0.1%. 
Unfortunately, our analysis is not sufficiently sensitive to determine the energy dependence 
of the parameters of the imaginary part. 

Figure [3] presents the predicted angular distributions calculated using this average local 
potential for each of the recently measured elastic scatte ring cross sections, together with 


the experimental data taken at University of Notre Dame [26|, to illustrate the overall good 


description. For completeness. Fig. [3] shows in the case of Ec.m. = 23.6 MeV and 26.0 MeV 
the local potential with additional volume contribution. The local hts using the parameters 
from Table [Hand the phase shift hts are not shown in Fig. [3] because they follow exactly the 
experimental data points (Xred ~ !)■ 


3.2. Comparison to predictions calculated with global a + nucleus optical potential 

In the framework of the 7 process network studies more than 2 x 10"^ reactions on 
about 2000 mostly unstable nuclei are taken into account. In many of those reactions alpha 
particles are present either in the entrance or in the exit channel. Experimentally, most of 
these reactions are inaccessible for cross section measurements, since short-lived unstable 
nuclei are involved in the entrance channel. Therefore, global a + nucleus optical potentials 
are used in the reaction network to predict the needed reaction rates. 

The variation of the potential parameters of the real part as a function of mass and 
energy is smooth and relatively well understood. On the contrary, the imaginary part of 
the optical potential is strongly energy-dependent especially at energies around the Coulomb 
barrier. As a result of these different energy dependencies of the parameterizations, different 
predictions for reaction cross sections will occur, in particular at very low energies far below 
the Coulomb barrier. 

In the following, we will compare our new experimental data to the predictions of well- 
known or recent open-access global potentials. 

(I.) In the middle of the 1960s, alpha elastic scattering experiments at Eq, = 24.7 MeV 
on nuclei ranging from O to U were performed and a consistent optical potential study 
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was carried out. Based on this study, a global a + nucleus potential was formed. It is 
a very simple, four-parameter Woods-Saxon potential with mass and energy-independent 
parameters referred to as the potential of McFadden and Satchler . This potential is the 
default for the H-F calculations of astrophysical reaction rates in the NON-SMOKER code 


191 146| . Despite its simplicity it provides a reasonable description of alpha scattering data 


and cross sections of a-induced reactions, at energies slightly above the Coulomb barrier. 
The H-F calculations based on this McFadden potential have a tendency to overestimate 


reaction cross sections at energies below the Coulomb barrier, e.g., [27|, 1^. Recently the 


depth of this potential was modihed by adding a Fermi-type function and this modihed 
potential provided excellent predictions for a-induced cross sections on and i 62 ,i 66 gj. 


221 24, 25 


(II.) The ’’standard” potential of M. Avrigeanu et al. ^ is an a -|- nucleus optical po¬ 


tential with Woods-Saxon parameterizations for the real and imaginary parts, with a total of 
nine independent parameters, three for the real potential, depth Vq, radius Rr and diffuse¬ 
ness qr, three to describe the imaginary volume potential Wy, Ry, o,v and the remaining 
three to describe the imaginary surface potential Ws, Rs, cts- All of these parameters are 
mass and energy dependent and the a -f nucleus optical potential was obtained by htting 
alpha particle elastic scattering and reaction cross sections on nuclei in the 45 < A < 209 
mass range, and for E < 50 MeV. Consequently, this potential provides a reasonable de¬ 
scription for both the scattering and reaction data. However, the extrapolation of such a 
many-parameter potential to unstable nuclei with extreme neutron-to-proton ratio N/Z may 
lead to additional uncertainties in the calculation of astrophysical reaction rates. 

In recent years, several a-induced cross sections on nuclei with mass number above A = 
140 0, lU, 0, in, lU Eil were measured and, furthermore, new experimental alpha elastic 
scattering angular distributions on even-odd nuclei became available j0| , thus improvements 
of this parametrization became possible. It was shown that the accuracy of the cross 
section predictions can be increased using an updated surface imaginary-potential depth 
Wr. However, for the ^°®Cd target under study in this work, the modihcations in the 
latest version are very minor and restricted to energies below about 10 MeV, i.e., the 
modihcations do not affect the anal ysis of the elastic scattering data. Therefore, we show 


the results using the potential from . 

(HI.) FinalW the most recent a -1- nucleus optical potential used in this work is the 
ATOMKI-Vl [2^ optical potential parameter set. It is based entirely on alpha elastic scat¬ 
tering angular distributions measured in the 89 < A < 144 mass region at energies above and 
below the Coulomb barrier. The real part V (r) consists of a double folding parametrization 
VA(r), based on the widely used DDM3Y interaction 29, 4l| and modihed by a strength 

parameter A derived from the energy-independent volume integral Jr. The Ws imaginary 
part is parametrized by a surface Woods-Saxon function with hxed Rs and as geometry and 
energy-dependent strength Ws and volume integral Jj. 

The ATOMKI-Vl potential has been constructed for low energies where a pure surface 
Woods-Saxon potential in the imaginary part is sufficient. Thus, it is not surprising that 
ATOMKI-Vl is not able to reproduce the angular distributions at the highest energies 
studied in this work. 
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Figure |3] shows the available experimental angular distributions in comparison with the 
predictions calculated using the above mentioned global optical potential parameter sets. 
It can be seen that each of the different a + nucleus optical potentials used in this work 
provides a slightly different description of the experimental data, however, none of them 
describes the entire experimental data set. For a strict comparison between the potentials, 
the values and total reaction cross sections cr^eac can be found in Table [H 


Table 2. xled analysis of the global optical potentials models considered in this work. Please note that the 
* indicates that the > 100 and ND refers to the Notre Dame measurements. 


Global Optical 


Elastic Scattering {Eq. 

M.) 

Potential Model 

15.6 MeV - Atomki 16.4 MeV - ND 

17.0 MeV - Atomki 

Avrigeanu et al. 2010 

4.3 

25.7 

2.7 

McFadden & Satchler 

35.1 

* 

41.9 

ATOMKI-Vl 

7.7 

1.8 

5.9 

Avg. local potential 

4.1 

1.7 

3.9 


17.0 MeV - ND 

18.3 MeV - ND 

18.9 MeV - Atomki 

Avrigeanu et al. 2010 

27.8 

22.1 

5.2 

McFadden & Satchler 

* 

* 

72.0 

ATOMKI-Vl 

1.6 

1.7 

1.3 

Avg. local potential 

2.2 

7.3 

4.9 


21.2 MeV - NDf 

23.6 MeV - NDf 

26.0 MeV - NDf 

Avrigeanu et al. 2010 

19.7 

20.7 


McFadden & Satchler 

* 

ik 

* 

ATOMKI-Vl 

27.2 

* 

* 

Avg local potential 

* 

★ 

-k 

t - Abs. renorm. 

1.016 

0.959 

0.921 


From the results of Table [H on the one hand, we found that the considered global a + 
nucleus optical potentials roughly describe the experimental data in the energy range from 
Ec.m. = 15.6 MeV to Ec.at. = 18.9 MeV. On the other hand, above Ec.m. = 18.9 MeV, all 
a + nucleus optical potentials fail to accurately describe the data, providing relatively large 
X^ed values. This asks for further improvement of global a-nucleus potentials but does not 
affect the following study of low-energy a-induced reaction cross sections. 


4. Comparison to CK-induced reaction cross sections on ^°®Cd 

Reaction cross sections of a-induced reactions for heavy nuclei can be calculated using 
the statistical model (StM) jsH. In particular, this model has been widely applied for the 
calculation of reaction cross sections and stellar reaction rates in nuclear astrophysics j^. 
In a schematic notation, the reaction cross section in the StM is proportional to 


a(a, X) ~ 


TgfiTx 


(3) 
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with the transmission coefficients Ti into the i-th open channel. The Ti are calcnlated from 
global optical potentials (particle channels) and from the 7 -ray strength fnnction for the 
photon channel. For details of the dehnition of Tj, see Ta^ refers to the entrance 

channel where the target nncleus is in its gronnd state. Note that this holds only for 
laboratory experiments; nnder stellar conditions, i.e., for the calcnlation of stellar reaction 
rates Na < ctv >, thermal excitations of the target have to be considered. 

It is typical for a-induced reactions on heavy nnclei that is mnch smaller than the 
other Ti. A simple qnalitative explanation is the high Conlomb barrier in the a channel. In 
the nentron channel, a Conlomb barrier is completely missing and, in the proton channel, the 
barrier is mnch lower. As a conseqnence, the cross section in the StM in Eq. ([3]) factorizes 
into a prodnction cross section of the componnd nnclens which is proportional to Tq, q, and 
a decay branching ratio bx = Tx/J2i'^i practically independent of T„. The prodnction 
cross section is thns entirely dehned by the chosen a + nnclens optical potential whereas 
the branching ratio bx depends practically not on the a + nnclens optical potential bnt 
on all the other ingredients of the StM (optical potentials for the other channels, q-ray 
strength fnnctions, level densities). Conseqnently, all cross sections of a-indnced reactions 
are sensitive to the a + nnclens optical potential, bnt each individnal (q!,p), or 

{a, 7 ) reaction has further and complicated sensitivities to the other ingredients. As the 
a + nucleus optical potential affects directly the production cross section, the sensitivity 
(as dehned e.g. in |^) is close to 1 for all (a,X) reactions at energies around or below the 
Coulomb barrier. 

The reaction Q-values for a-induced reactions on ^°®Cd are —10.15MeV for the (a,n), 
—5.51MeV for the (a,p), and -1-1.13MeV for the ( 0 , 7 ) reactions. Experimental data for 
these reactions are available at low energies from about 7 MeV to 13 MeV . Following 
the above argumentation, it can be expected that the [a, n) cross section is dominating at 
higher energies down to close above the {a,n) threshold. Below 10 MeV, only the (a,p) 
and ( 17 , 7 ) channels are open, and it turns out that the {a,p) channel becomes very weak 
towards lower energies. 

The focus of the present study is the a T nucleus optical potential for ^°®Cd which 
determines the compound formation cross section. In Fig. |^, we show the formation cross 
section taken as the sum over the dominating ( 0 , 7 ), {ct,n), and {a,p) channels for the 
potentials under study in this work. Obviously, the calculated compound formation cross 
section should be an upper envelope of the experimental ( 0 , 7 ), (a,n), and {a,p) data. At 
higher energies above about 13 MeV, all potentials lead to very similar cross sections. This 
is a general Ending which has been explained from the properties of reflexion coefficients 
Pl, see ( 53 . 55l |. However, at lower energies, significant differences can be seen. The local 
potential which has been adjusted to the elastic scattering angular distibutions overestimates 
the experimental data at energies below about 10 MeV. Unfortunately, it was not possible to 
determine the energy dependence of the imaginary part from elastic scattering and, therefore, 
the local potential behaves similar to the widely used potential by McFadden and Satchler 


45l | which has also been determined from elastic scattering at higher energies and also has an 


energy-independent imaginary part. The reaction cross sections from the McFadden/Satcher 
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potential have already been shown in 

The new ATOMKI-Vl potential uses a very similar parametrization as the local poten¬ 
tial, but it has an energy-dependent imaginary part in addition. The decreasing imaginary 
strength towards lower energies leads to smaller cross sections, but the calculation still over¬ 


estimates the experimental data. On the contrary, the potential by Avrigeanu et al. [48 
underestimates the data. The latest revision of this potential ji^ has a slightly stronger 
imaginary part at very low energies which leads to slightly larger cross sections by about 
10% below lOMeV. 

Let us now turn to the individual ( 0 , 7 ), (Q!,n), and (a,p) reactions. As already stated 
above, these cross sections depend not only on the a + nucleus optical potential with 
sensitivity ^ ~ 1 as dehned in j^, but also on the further ingredients of the StM with 
varying energy-dependent sensitivities S. This complicates the comparison of the results 
of the StM with the experimental data. The following calculations are the results from 
SMARAGD code with its default parameters. Following the scope of this study, only 
the a + nucleus optical potential has been varied. 

Fig. O shows the astrophysical S factors of the ^°®Cd(a, 7 )^^°Sn reaction which domi¬ 
nates below 10 MeV. Here the (a, 7 ) cross section is essentially dehned by the a -|- nucleus 
optical potential. Similar to the compound formation cross section in Fig. |H the local po¬ 
tential overestimates the low-energy data signihcantly, the ATOMKI-Vl potential slightly 
overestimates, and the potential by Avrigeanu et al. slightly underestimates. At higher en¬ 
ergies above the neutron threshold, the cross sections become < 7 ( 0 !, 7 ) ~ TaflT.y/ with 

Ti ~ Tn and, thus, no direct conclusion on the a + nucleus optical potential is possible 
here. As a word of caution, it must be mentioned that the behavior of the {a, 7 ) cross section 
at low energies may also be affected by direct reactions, in particular. Coulomb excitation 


of the hrst 2+ state in ^°®Cd 57 


Fig. O shows the same comparison for the ^°®Cd(a, n)^‘’®Sn reaction. At energies suffi¬ 
ciently above the threshold, cross sections become the a{a,n) Ta,oTn/ Y.i Ti ~ Tafi for 
Ti ~ Tn- In principle, it would be possible to derive properties of the a - 1 - nucleus optical 
potential from these data; however, the experimental data cover only a narrow energy in¬ 
terval of about 2 MeV from 10 MeV to 12 MeV, and the error bars are not sufficiently small 
to allow clear conclusions. 

Fig. [7] shows the astrophysical S factors of the ^°®Cd(a,p)^‘^®In reaction. Similar to the 
(a,n) data, the energy range is also very limited (9MeV - 10MeV). As the cross section 
becomes a{a,p) ~ Ta^Tp/ J2i Ti ~ TapTp/(Tp + T.y) in the energy range of the experimental 
data, it is not possible to directly conclude on the a + nucleus optical potential from these 
data. 


^In astrophysical investigations it is common to quote the astrophysical S factor . The cross section, 
cr(E) and the astrophysical S factor, «S'(E) at c.m. energy E are related by: 

S{E) = Ea{E)exp{2Trri), (4) 

with 1 ] being the Sommerfeld parameter. For better visibility, in the figures instead of the cross sections, 
the astrophysical S factors will be used throughout the present manuscript. 
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For completeness it should be mentioned here that the study of cross section ratios like, 
e.g. cr(a,p)/cr(a, 7 ) ~ {TaflTp/ Y,^Ti)/[TaflT^/ = Tp/T^ allows to constrain the other 

ingredients of the StM except the a + nucleus optical potential. For an application of this 
idea, see e.g. j^. As we have seen (and has been stated in our earlier study of the reaction 
data [^), it is not straightforward to draw conclusions on the a + nucleus optical potential 
from the cross sections of the ( 0 , 7 ), (a,n), and (a,p) reactions. 


5. Summary and conclusions 

We have measured three angular distributions of ^°®Cd(a, Q!)^°®Cd elastic scattering at 
energies around the Coulomb barrier; and we have analyzed these data and data from the 
literature within the optical model. Reasonable descriptions can be obtained for all data 
using a real folding potential with two adjustable parameters (strength A and width w) and 
an imaginary part of Woods-Saxon type. Similar to previous results jlH, at low energies, a 
surface Woods-Saxon type parametrization is sufficient, whereas at higher energies a com¬ 
bination of volume and surface imaginary part improves the obtained signihcantly. A 
local potential is derived from the average parameters. Unfortunately, it was not possible 
to determine the energy dependence of the imaginary part at low energies. The new exper¬ 


imental data are also described very reasonably by the recent global potentials [28], 1^, [49 
whereas the old McFadden/Satchler potential |4^ has a trend to overestimate the elastic 
data. 

The locally adjusted potential and the global potentials by Avrigeanu et al. |48| and 
ATOMKI-Vl 1281 were used to calculate the cross sections of a-induced reactions on ^°®Cd. 


It is shown that these cross sections factorize into a compound production cross section which 
depends only on the chosen a + nucleus optical potential and into a decay branching which 
is practically independent of the a + nucleus optical potential but shows a complicated 
dependence on the other ingredients of the StM. Therefore, clear conclusions on the a 
+ nucleus optical potential from a-induced reaction data are difficult in the case of ^°®Cd. 
Nevertheless, some conclusions can be drawn. It turns out that the locally adjusted potential 
overestimates the reaction data at low energies. This result is similar to the one obtained 
with the widely used McFadden/Satchler potential and can be attributed to the missing 
energy dependence of the imaginary part. The predictions from the new global potentials j28|, 
^1 are much closer to the experimental data but still far from an excellent description 
with Xred ~ 1- The ATOMKI-Vl potential slightly overestimates the reaction data at 
low energies whereas the potential by Avrigeanu et al. underestimates the data. As similar 
results have been found for ^'^°Ce and ^^^Pr, the best prediction of a-induced reaction cross 
sections in the mass range 100 < A < 150 could be taken from the geometric mean of these 
recent global potentials. 
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0c.m, (deg) 0^,^, (deg) 


Figure 3. (Colour online) Available experimental ^°®Ccl(a, a)^°®Cd elastic scattering cross sections — 
measured at Atomki and at University of Notre Dame [2^ —, normalized to the Rutherford cross section 
compared to theoretical predictions calculated using global a + nucleus optical potentials. For details, see 
text. 
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^°®Cd(a,X) 



Figure 4. (Colour online) Total reaction cross section for a + from different potentials, converted 

to the astrophysical S-factor, and compared compared to the experimental (a, 7 ), {oi,p), and (a,n) cross 
sections [l^- Further discussion, see text. 


^°®Cd{a,y)"°Sn 



Figure 5. (Colour online) Astrophysical S-factor for the ^°®Cd(a, 7 )^^°Sn reaction, experimental data from 
[27| . The figure contains the predictions from global a + nucleus potentials 28|, and the average local 
potential. 
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Figure 6. (Colour online) Astrophysical S-factor for the ^°®Cd(a, n)^°®Sn reaction, experimental data from 
[27| . The figure contains the predictions from global a + nucleus potentials 28|, and the average local 
potential. 





Figure 7. (Colour online) Astrophysical S-factor for the ^^°Cd(a,p)^^^In reaction, experimental data from 
[271 1. The figure contains the predictions from global a 
potential. 


nucleus potentials [28], |48| and the average local 
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